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Antiferromagnetic crystalline topological insulators 
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The gapless surface Dirac cone of time reversal invariant topological insulators is protected by 
time reversal symmetry due to the Kramers' theorem. Spin degree of freedom is usually required 
since Kramers' theorem only guarantees double degeneracy for spinful fermions, but not for spinless 
fermions. In this paper, we present an antiferromagnetic spinless model, which breaks time reversal 
symmetry. Similar to time reversal invariant topological insulators, this model possesses a topolog- 
ically non-trivial phase with a single surface Dirac cone, which is protected by the combination of 
time reversal and translation operation. Our results show that in magnetic crystals, a single Dirac 
cone can exist on the surface even without any spin degree of freedom and spin-orbit coupling. 



Time reversal (TR) invariant topological insulators 
(TIs) [l|-y| possess insulating bulk states, and metallic 
surface states consisted of an odd number of Dirac cones, 
e. g. a single Dirac cone in the prototype TIs, Bi2Se3 
family of materials [5|-|7| . The gapless nature of topolog- 
ical surface states is protected by TR symmetry accord- 
ing to the Kramers' theorem |8l4ll|. For a TR invariant 
system, TR operation 6 satisfies 6^ = — 1 for spinful 
fermions and O^ = 1 for spinless fermions or bosons. Ac- 
cording to the Kramers' theorem, when 0^ = — 1 there is 
a double degeneracy for each energy level, which protects 
the gapless nature of surface Dirac cones in TR invari- 
ant TIs. Therefore, spin degree of freedom, as well as 
spin-orbit coupling, is required for TR invariant TIs. 

The idea of topological phases due to the symmetry 
protected degeneracy is generalized to other types of 
symmetry, such as crystalline symmetry 12M2Q] , particle- 
hole symmetry 2ll423j , etc. In this work, we present a 
spinless antiferromagnetic model with topologically non- 
trivial phases, which breaks TR symmetry, but possesses 
surface states of a single Dirac cone, similar to that of 
TR invariant TIs. The combined symmetry of transla- 
tion operation and TR operation plays the role of TR 
symmetry in TR invariant TIs and protects the gapless 
nature of Dirac cones, which is similar to the previous 
work of R. Mong, et at [24]. Different from that work, 
our model is spinless so that the TR operation by itself 
satisfies O^ = 1 and cannot protect any degeneracy. The 
translation operation plays an essential role in inducing 
non-trivial topological phases, and in this sense, we dub it 
as "antiferromagnetic crystalline topological insulators" . 
Therefore, our model indicates that spin and spin-orbit 
coupling are not necessary for topological phases with a 
single surface Dirac cone in magnetic structures. 

We study a model for a layered antiferromagnetic 
structure with each layer consisted of a square lattice 
with magnetic moments on each lattice site perpendicular 
to the xy plane, as shown in Fig. [T]^b). We consider the 
case that magnetic moments in one layer are polarized 
along the same direction, so each layer is ferromagnetic. 
But the magnetic moments between two adjacent layers 
are polarized in the opposite directions, so the whole sys- 



tem has an antiferromagnetic structure with two atoms 
in one unit cell, denoted as A and B, as shown in Fig. 
[H (a) and (b). In one layer, there are three orbitals |s), 
\Px) and \py) on each lattice site. The \px) and \py) or- 
bitals carry the angular momentum 1 and couple to mag- 
netic moments through Zeeman type of coupling. In one 
layer, the lattice vectors are denoted as ai = (a, 0, 0) and 
^2 = (0, a, 0), while layers are stacked along the direction 
as = 2c = (2ci, 2c2, 2C3), where the vector c points from 
the A atom to the adjacent B atom, as shown in Fig. [1] 
(a). The Hamiltonian is given by 

H = H+^H-^ He, (1) 

77 / V nm anr] P^^^ ' / v " anr] otnr] 

{n'm)in,a(3 ft, a 

+ Yl ^^1 (-*4p.^cfip„'; + *4p.^cHj,^^) (2) 



^c = Yl (^Hi4a+c/3m- + h.c.) 

{nfh) ciOcl3 



(3) 



where 77 = + is for A layers and 77 = — for B lay- 
ers, n = {ux^riy^nz), rh = {rrix^my^mz) denote lattice 
sites and a, /3 = s,Px-,Vy denote orbitals. {nm)in rep- 
resents the nearest neighbor hopping in xy plane with 
the hopping parameters t^^ while {nm)c indicates the 
nearest neighbor hopping between A and B atoms along 
c axis with the parameters r^^. We take into ac- 
count the (J bond for s orbitals, the a and tt bonds 
for p orbitals, and the a bonds between s and p or- 
bitals. For the in-plane hopping, t^^ are parameter- 



ized as tf^^, = t%%^^ = tf^, ^ = t%-^, ^ = Usa, 
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spa , where Cx and Cy denote the 



n-\-ey,n 

vectors to the nearest neighbor site along x and y direc- 
tions, respectively. For the hopping between two layers, 
since c is not along z direction, we need to decompose 
p orbitals into components along c axis and perpendicu- 
lar to c. Consequently, we obtain rl^^^^ = r^^^^^ = Vga, 



n,n-\-c 



n+c,n ~ ^pcr^l^^pTvy-i^ ^l), ' fi,n+c ~ ' n+c,n 
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^A2, where Ai 



(A2 = f) is 



— £1. 

the angle between c and x (y) axis. Mi term is the Zee- 
man type of couphng between p orbitals and magnetic 



J 



moments. In the momentum space, the Hamiltonian is 
given by 



Hr, 



^^1^ I 2msp^sin(27r/ci) 
2iUspa sin(27r/c2) 



-2iuspa sin(27r/ci) —2iuspa sin(27r/c2) 
-ir]Mi 



Epx\J^) 
ir]Mi 



^r 
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where 



Es{k) = 2usa{cos{27rki) -\- cos(27r/c2)) + e^ 
Epx{k) = 2(i^pcr cos(27r/ci) + 2x^7^ cos(27r/c2)) + « 



Eny\k) 



^py\ 



2(iip7r cos(27r/ci) + iipcr cos(27r/c2)) + e^. 



Here 






(ct^(^),4^^(^),4^^(^)), 



and k 



(6) 



Si=i 2 3 ^«^* with the reciprocal lattice vectors hi sat- 
isfying 6i • dj = 27T5ij. The basis wavefunction for the 
above Hamiltonian is given by 

1 X- 



\ar],k) 



k-fj 
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^0a(r- 



' nr] ) 
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where A^ is the normalization factor, fnr] = Rn + ^77 with 
the lattice vector Rn and the position hr^ for the atom A 
or B in one unit cell, and (j)cx denotes the basis wavefunc- 
tion for the orbital a {a = s^Px^Py). 

When the coupling Mi is large enough, each layer 
will be driven into the quantum anomalous Hall phase 
with a non-zero integer Hall conductance. To see this, 
we may consider one layer with the effective Hamilto- 
nian (|4|) and change the basis from \px) and \py) to 
|P+) = -^i\Px) -^^Py)) and |p_) = -^HPx) - i\Py))- 
The effective Hamiltonian (|4]) is rewritten as 



Hf] = 



— V2iu, 



tUs 



Es{k) 
sp^(sin(27r/ci) 
pa{sm{27Tki) 4 



>/2mspcr(sin(27rA:i) + i sin(27r/c2)) —V2iuspa{s^^{27rki) — i sin(27r/c2)) 



- i sm{27Tk2)) \{Ep^(k) + Epy{k)) + iqM^ 
zsin(2^/c2)) -\{Ep^{k) - Epy(k)) 



-\{Ep,{k) - Epy{k)) 

\{EpS) ^ Epy{k)) - r]Mi 



If the parameters of this model satisfy the condition 
\Es{k)-Ep{k)-Mi\ < 2|Mi| {Ep = \{Ep,^Epy)), s or- 
bital and one of the p orbitals with the energy Ep{k)-\-Mi 
are close in energy for both A and B layers, while the 
other p orbital are well separated. Let us assume the 
Fermi energy is between s orbital and the p orbital with 
the energy Ep{k) + Mi, and then the low energy physics 
is described by the two-band effective Hamiltonian 



H, 



effjT] 



s{k)- 



E 

^x,y,z 



djCFi 



(9) 



where <j is the Pauli matrix denoting the basis of s and 
p orbital. For A (B) layer r] = ±, e = \{Es(k) + Ep{k) + 



Ml), dx 
and dz = 



-V^iispcrSin(27r/c2), dy = ^v^TXspcr sin(27r/ci) 
Ep{k) - Ml). This model has been 
and it possesses a non-zero quan- 
tized Hall conductance when the configuration of the vec- 
tor d = ^{dx^ dy^dz) {d = \/^i df) in the whole Brillouin 



\{Es{k) 
well studied in Ref. IZ 



I 

zone has a non-zero winding number. For example, we 
can consider the parameter regime when this model has a 
direct band gap around F point (/c = 0). Expanding the 
Hamiltonian ([9j) around F point, we find the low energy 
physics is described by Dirac Hamiltonian with linear dis- 
persion. For Dirac Hamiltonian, a transition happens at 
£^s(0) = Ep({)) + Ml when the band gap closes. Since 
the system is a normal insulator for Mi = 0, both lay- 
ers are driven into the quantum anomalous Hall phase 
with chiral edge states when |Mi| > \Es{^) - Ep({))\ , 
as depicted in Fig [1] (b). Furthermore, because of an- 
tiferromagnetism, the Hall conductances, as well as the 
velocities of the corresponding chiral edge states, are op- 
posite for two adjacent layers. Therefore, one may expect 
the tunneling between two layers will destroy the chiral 
edge states and lead to a gap opening for surface states. 
To test this idea, we perform numerical calculation of the 
energy dispersion for both the bulk and slab configura- 
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FIG. 1. The schematic plot of the lattice structure, (a) For 
each layer, it is the square lattice with three orbit als, s, px 
andpy orbitals. The A layer is shifted by the vector c relative 
to the B layer, (b) The lattice is formed by the stacking of 
A and B layers and each layers can be viewed as the QAH 
layers. 
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tion, as shown in Fig. [2] (a) and (b). Here we take the 
vector c = (0,0,C3) and consider the side surface normal 
to y direction for the slab calculation. From Fig. [2] (a), 
there is indeed an energy splitting for the bulk dispersion 
at the V — X — Y — M plane (see the Brillouin zone in Fig 
[2](c)), but the energy levels at the Z — U — R—T plane are 
still degenerate. Correspondingly, surface states open a 
gap at r point, but remain gapless at Z point and form a 
single Dirac cone in the slab configuration. Here f and Z 
are the projection of F and Z in the two-dimensional sur- 
face Brillouin zone, respectively. In the following, we will 
show that the gapless point of surface states at Z point 
is protected by the symmetry operation Tc9, where Tc is 
the translation operation from A layer to B layer along 
the vector c and O is TR. 

In TR invariant TIs, TR symmetry O plays an essen- 
tial role. For spinful fermions, = iciyK where K is the 
complex conjugate and ay is the Pauli matrix for spin, 
so ©^ = —1, which guarantees the double degeneracy 
for all the states at TR invariant momenta in Brillouin 
zone following Kramers' theorem. This double degen- 
eracy protects the gapless point of topological surface 
states. In constrast, for spinless fermions, = K and 
6^ = 1, so TR can not protect any degeneracy by iteself. 
The present model has ferromagnetism within one layer 
and antiferromagnetism between two adjacent layers, so 
the Hamiltonian has a combined symmetry TcQ. To un- 
derstand the role of Tc 6, we may calculate the behavior 
of (Tc6)2 on the basis (O and find 






(10) 



-T] means that A and B layers are inter- 
depends on the value 



ai^Tvks 



where 77 

changed. Thus, (T^Bj^ 
of ks. For the plane ks = 0, (T^Oj^ = 1, but for the 
plane fe = ^, (Tc6)^ = —1, which explains the dou- 
bly degenerate states at TR invariant momenta when 



FIG. 2. (a) Bulk band dispersion of the present model with 
the parameters given by Usa — —0.2, Uspa = 0.2, Upa = 0.2, 

UpTT = 0.2, Vsa = 0.05, Vspa = 0.3, Vpa = "0.1, Vp^ = 0.1, 

Cs = 0, e = -5, Ml = 4.5, a = 1 and c = (0,0,2). (b) The 
energy dispersion for a slab of the present model with the 
same parameters as above. A Dirac cone appears at the Z 
momentum, (c) The Brillouin zone of the lattice and the TR 
invariant momenta in the Brillouin zone. 
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In Fig [2] (a), we find there is a double de- 
generacy for each momentum in the ks = ^ plane and 
this is because we take c= (0, 0, cs) along z axis and the 
Hamiltonian has additional two-fold rotation symmetry 
C2 around z axis. On the basis \ar]^k)^ the C2 opera- 
tion is written as C2\o:r],i^,ks) = i^^lar^, — /^, /cs), where 
^ = (/ci,/c2) and Ua is 1 for s orbitals and -1 for p or- 
bitals. Therefore, for any eigen state |^rn(^)) in the 
tv = ^ plane, where m denotes different eigen states, 
C2TcO|^rn(^)) is also an eigen state with the same eigen 
energy at the same momentum. Moreover, when /C3 = ^, 



TcO|a±, /^) 



}b3-c-\-ib3-hfj 



\aT,-i^) 



^z|a^, — ^), so 



one can see that the operation TcO at the ks = ^ plane is 
quite similar to TR operation for spinful fermions. Con- 
sequently, one can show that (^rn(/^)|C'2Tc9|^rn(/^)) = 

(ec2Te6^^(^)|e^^(^)) = -{^mmc2Tce\^mm = 

0, so the eigen states C2Tc6|^rn(/^)) and |^rn(/^)) are or- 
thogonal and there is a double degeneracy for each mo- 
mentum /^ at /c3 = I plane. 

When ci,2 in c is finite, C2 rotation symmetry is bro- 
ken, so the degeneracy at /C3 = ^ plane is removed, ex- 
cept four momenta Z^ U^ R and T, as shown in Fig [3] 
(a). These momenta are invariant under TcB, so TcQ 
by itself protects the degeneracy. One can directly show 
(^^(A)|Tee|^^(A)) = -(^^(A)|Tc6|^^(A)) = with 
A = Z, /7, R, T, so l^m(A)) and Tc6|^m(A)) give two de- 
generate eigen states. From Fig. [3](b), one finds that the 
gapless Dirac point of surface states at Z also remains. 



The Z2 topological invariant can be defined in the 
ks = ^ plane due to the degeneracy at Z, U^R^T. Since 
we have only four Tc 0-invariant momenta, only one Z2 
number is allowed, and given bv|ll| 



(-i)-^ = UA 



Si = 



^det[w{Ai)] 
Pf[w{Ai)] 



(11) 



where u is Z2 topological number and A^ is taken to 
be Z, /7, R, T. The matrix w is defined as Wmn{f^) = 
{^rn(— /^)|^c0|^n(/^)) with m and n taken from all the 
occupied bands. det[w{Ai)] and Pf[w{Ai)] are the de- 
terminant and pfaffian of w at A^. Alternatively, one 
can also define Z2 invariant with the Berry gauge poten- 



tial A,,^ = Em (^ 



xy 



d A d_ A 

dk^ ^y dky ^x 



I'^rn) and the Berry curvature 
where m is for all the occupied 



bands. The Z2 number is given by [26] 



27r 



dn • A— I d^nFxy 

0\l\ /<2 ^ ^1/2 



mod2 (12) 



where 1^1/2 takes half of the /^s = | plane. The expression 
(p!2]) is much easier to evaluate numerically and with the 
numerical method introduced in Ref. [27], we obtain v = 
1 for the parameters given in the caption of Fig. [2] and 

13 

For TR invariant TIs, when there is an inversion sym- 
metry, Z2 invariant can be easily determined by the par- 
ity eigenvalue of the occupied bands [11]. However, this 
simplification cannot be applied here because inversion P 
does not commutate with T^G, PT^B = T_c6P ^ T^BP. 
It turns out that C2 symmetry, which commutates with 
TcB, plays the role of inversion symmetry in TR invari- 
ant TIs. If there is C2 symmetry in our model, which 
requires the vector c along z direction, Z2 invariant can 
be calculated by the eigenvalues of C2 rotation of all the 
occupied bands, with the expression [11 1 



(-ir 



n,(5,: 



(5i = n-=i6m(Ai) (13) 



where C2m(Ai) = ±1 is the C2 eigenvalue of the 2mth 
occupied band. The total number of the occupied bands 
is 2N . Since the (2vn — l)th and 27Tith bands are related 
to each other by TcB, these two states share the same C2 
eigenvalue, ^2m(Ai) = ^2m-i(Ai). For our present model, 
s orbitals are even under C2 rotation while p orbitals 
are odd. Therefore, topological phase transition happens 
when there is a band inversion between s orbital and 
one of the p orbitals at high symmetry momenta. This 
is consistent with the above analysis based on the two- 
band effective model. With the parameters given in the 
caption of Fig [2l one find i^ = 1 for our model, which 
confirms the topological nature of gapless surface states. 
The model presented here is quite similar to the anti- 
ferromagnetic model in Ref. [2^, but there is one essen- 
tial difference. Ref. [24] discusses the model of spinful 
fermions with the TR operator B = icFyK^ so B^ = —1 by 
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FIG. 3. (a) Bulk band dispersion with the lattice vectors 
a—\ and c— (0.1, 0.15, 2). All the other parameters are the 
same as above. The degeneracies at Z, U, R, T are preserved 
due to the TcG symmetry, (b) Energy dispersion for a slab of 
the present model. 



itself. Our model is for spinless fermions, so B^ = 1 and 
the degeneracy originates from the combination of time 
reversal and translation operation. Consequently, in Ref. 
0, (Tc^f = -1 for ks = plane and (TcB)^ = 1 for 
ks = ^ plane, yielding the Z2 topological invariant only 
well-defined at /C3 = plane. It is exactly opposite in 
our case. In Ref. [24], for the Z2 topological invariant at 
the ks = plane, there is no difference between inver- 
sion P and two- fold rotation C2 along z axis. But in our 
case, these two operations are different for the ks = ^ 
plane, and only C2 rotation can be used to define Z2 
invariant. Since Z2 invariants are defined at different 
momentum planes for two cases, the corresponding sur- 
face Dirac cones also appear at different high symmetry 
momenta. Our model shows that for magnetic systems, 
topological phases with a single surface Dirac cone can 
exist in a spinless fermion model. 

In conclusion, we study a spinless model in an antifer- 
romagnetic structure, which has the topologically non- 
trivial phase with a single surface Dirac cone. The bulk 
Z2 invariants are evaluated and compared with the direct 
numerical calculation of surface states in a slab configu- 
ration. Our model indicates the possibility to generalize 
the concept of TIs to magnetic structures [28| and other 
boson systems. 

We would like to thank C. Fang, X. Liu, Y.H. Wu for 
useful discussions. 
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